Effective Edwards-Wilkinson equation for single-file diffusion 
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In this work, we present an effective discrete Edwards- Willcinson equation aimed to describe tlie 
single-file diffusion process. Tlie fcey physical properties of the system are captured defining an 
effective elasticity, which is proportional to the single particle diffusion coefficient and to the inverse 
squared mean separation between particles. The effective equation gives a description of single-file 
diffusion using the global roughness of the system of particles, which presents three characteristic 
regimes, namely normal diffusion, subdiffusion and saturation, separated by two crossover times. 
We show how these regimes scale with the parameters of the original system. Additional repulsive 
interaction terms are also considered and we analyze how the crossover times depend on the intensity 
of the additional terms. Finally, we show that the roughness distribution can be well characterized 
by the Edwards- Wilkinson universal form for the different single-file diffusion processes studied here. 

PACS numbers: 05.60.-k,02.50.Ey,05.40.Fb 



o 

CD 



I 

o 
o 



> 
00 

in 
o 

\6 
o 
o 



X 



I. INTRODUCTION 

Single file diffusion (SFD) is one of the simplest and 
non-trivial problems in statistical physics In this 

problem a one-dimensional collection of particles diffuses 
with hard-core repulsion as the basic constraint. In one 
dimension, hard-core repulsion implies the non-passing 
rule, i.e. a given particle can not cross over its neigh- 
bor particles. This restricts the phase space generating 
an inherent entropic interaction. Despite its simplic- 
ity the SFD problem still attracts much attention due 
mostly to the possibility of applying different analytical 
formalisms [^-Q- 

A direct experimental observation of SFD processes 
has become possible during the last decade. The im- 
provement of experimental techniques has permitted 
studying systems ranging from interacting colloidal parti- 
cles on channels or rings at the microscopic scale 
to charged metallic balls on rings at the macroscopic 
scale [1^. Other interesting problems in which single- 
file diffusion is important are the longitudinal relaxation 
of a single polymer fisl IT^ or the transverse fluctuation 
of vecinal surfaces [15| . 

The simple hard-core interaction gives rise to non- 
trivial diffusion properties. If the diffusion of a tagged 
particle is followed, it first presents a normal diffusion 
regime, i.e. A(t) = {[r{t) - {r){t)Y) = 2Dt, where r{t) 
is the instantaneous position of the particle, (...) is an 
ensemble average, and D is the diffusion coefficient. This 
relation is valid up to a characteristic time, beyond which, 
in average, particles start colliding with their neighbors. 
Then, as a consequence of the interaction between par- 
ticles and the one-dimensional character of the problem, 
the tagged particle shows a subdiffusive behavior of the 
form A{t) = 2Fii/2, where F is a characteristic coef- 
ficient. This subdiffusive regime last for ever in infinite 
systems, but it crosses over to saturation or normal diffu- 
sion in finite systems, depending on the applied boundary 



conditions [1] . The subdiffusive regime of a single tagged 
particle is the signature of the interacting system. On 
the contrary, the diffusion of the center of mass is al- 
ways normal, i.e. it increases linearly with time, without 
showing any signature of the interacting system. 



Different analytical approaches to the SFD problem 
have been recently proposed. Among others, we can men- 
tion a Langevin formulation [8], the calculation of the 
probability density function for the position of a tagged 
particle using a Bethe ansatz [6, 7], an approach based on 
fractional Brownian motion with a variable Hurst expo- 
nent [l6j , or a formulation for the tagged particle motion 
using classical reflection and transmission coefficients @ . 
In the present work, we propose a different and comple- 
mentary approach based on an effective description of the 
SFD problem in terms of a discrete Edwards- Wilkinson 
equation (DEWE). The two parameters of the DEWE 
are the noise intensity and the elasticity of the system. 
In order to present an effective equation, effective val- 
ues for these parameters are given, relating them to the 
SFD parameters, i.e. to the mean particle separation 
and to the single particle diffusion coefficient. Then, we 
describe the SFD process using the global roughness of 
the system of particles, which gives information on how 
the particles diffusion becomes correlated with time. It 
is shown that the description using the DEWE success- 
fully accounts for the case of a pure hard-core interac- 
tion. In order to extend our analysis, we consider two 
further cases of interacting particles, one with a purely 
harmonic repulsion term, and the other with an interac- 
tion term proportional to the reciprocal of the squared 
distance between two consecutive particles. Finally we 
discuss some aspects of the universal properties of the 
distribution function of the global roughness. 



II. MODEL SYSTEM 

Let us first briefly describe tlie model system under 
study. We consider a set of N particles of unitary size 
in a discrete one-dimensional system. Tlie particles are 
in average a distance L apart and thus the size of the 
system is LN and the particle density is p = 1/L, as 
depicted in Fig. [T] For the sake of simplicity, we use 
periodic boundary conditions. The j label is used for 
the particles, so the instantaneous position for the j-th 
particle at time t is rj(t). Instead of this instantaneous 
particle position, let us consider the displacement field 
Uj(t) = rj{t)—jL, which measures the displacement with 
respect to an ideal ordered lattice rj = jL. 

We shall consider the motion of the particles in a dis- 
crete medium. The dynamics of each particle consists 
of a random displacement to one of its two neighbor- 
ing sites, with equal transition probabilities r„^„+i — 
r^.u-i = r = 1/2, only if the new site is not occupied 
by one of the neighboring particles. This means that the 
hard-core repulsion restricts the position of the j parti- 
cle to rj-i{t) < rj{t) < rj+i{t) for all times. This rep- 
resents indeed an entropic repulsion between particles. 
Note that, expressed in terms of the displacement field, 
the non-passing rule may be written as 

Uj-i{t) - L < Uj{t) < Uj+i{t) + L. (1) 

Within this simple model, the interacting particles fol- 
low a SFD process. Thence the single particle diffusion 
shows three characteristic regimes: a single particle nor- 
mal diffusion regime A = {[rj{t) — {rj{t))]'^) ^ t, an in- 
teracting particles subdiffusive regime A ~ i^/^, and a 
center-of-mass normal diffusion regime A ^ t. The sub- 
diffusive regime is the key signature of the interaction 
between particles. 

As described here, this simple model with hard-core 
repulsion and constant diffusion coefficient F does not 
depend on temperature. One can think of other inter- 
action terms which allow for a temperature dependence. 
For instance, by including a simple harmonic repulsion 
between particles given by the Hamiltonian 

Periodic boundary conditions are implicitly assumed in 
the sum. We shall analyze the model containing the hard- 
core repulsion plus the harmonic contribution as the sim- 
plest model with an energetic (non-entropic) interaction 
term. We shall also consider another energetic interac- 
tion term where the interaction is proportional to 
being r the distance between two consecutive particles. 
We are interested in this interaction form since it is rel- 
evant for the vecinal surface problem [l5l |. In such case, 
the Hamiltonian reads 

Af-l 

Ha^AJ2 K+1 (*) - % W + L]''- (3) 

j=0 
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FIG. 1: (Color online) Schematic representation of the model 
single-file difTusion system. Each particle of unitary size 
diffuses in a discrete one-dimensional system with periodic 
boundary conditions. The mean separation between particles 
is L. Each particle can move with equal transition probability 
r = 1/2 to one of its next nearest neighbor sites if they are 
not occupied by other particle. 



In order to compare with the effective DEWE descrip- 
tion of the process we present in the following results from 
direct numerical simulations of the present model. A 
Monte Carlo scheme can be implemented to perform nu- 
merical simulations of the model described so far. For the 
single-file diffusion mechanism a random particle is se- 
lected and one of its two neighbors randomly chosen with 
equal probability, and then the particle moves only if the 
new site is empty, which implements the non-passing rule 
Eq. ([TJ. When an energetic contribution is considered, 
it is taken into account in the dynamics of the particles 
with a Metropolis rule. The new site for a given particle 
is randomly selected as before, then the energy difference 
AH is considered and the probability P = exp(— AiJ/T) 
is computed. The new particle moves to the new site 
if P > ^, where ^ is a random number uniformly dis- 
tributed in (0, 1]. This Monte Carlo scheme will be used 
to compare with the effective description of the model in 
terms of the DEWE, which is presented in the following 
section. 



III. DISCRETE EDWARDS- WILKINSON 
EQUATION 

In this section we present the effective DEWE which 
allows for a description of the main physical properties 
of the SFD problem. We shall first present the well- 
known general solution of the equation and we will then 
discuss how its parameters can be effectively related to 
the parameters of the interacting particle system. Notice 
that although we use here a continuous-time Edwards- 
Wilkinson equation for the sake of simplicity, the results 
we will present here can be directly translated to the 
discrete-time version. 



A. General results 

The Edwards- Wilkinson equation is characterized 
by the amplitude of the thermal noise D and the elastic 
constant v. In a discrete form it can be written as 

dtUj{t) = [uj+i{t) + ~2uj{t)] + T]j(t), (4) 
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FIG. 2: Roughness of the DEWE for parameters = 128, 
1/ — 1 and D — 1/2. Dashed lines correspond to the three 
asymptotic regimes, which are separated by the crossover 
points yi = {ti,Wi) and 1/2 ~ {t2,Wi), as shown in the 
figure. 



with j = 0, — 1 and periodic boundary conditions. 
The noise term has a zero mean and the correlations 
{r]jit)r]j,{t')) = 2DSjj,S{t-t') have intensity D. The av- 
erage over the white noise is denoted by angular brackets. 
Notice that we write the DEWE directly in terms of the 
displacement field variable, thus elastically connecting a 
system of N interacting beads. We consider the center 
of mass of the displacement field as 



N 

which serves to define the Fourier transform 



J=0 



(5) 



(6) 



with fc„ — 2TTn/N. Then the DEWE is transformed to 
a set of independent equations for the time evolution of 
the Fourier modes, 



with 



dtCn{t) ^ -VnCn{t) + Vnit), 



i>n = sin^(7rn/A^). 



(7) 



(8) 



The sinusoidal character of this damping constant is asso- 
ciated to the discreteness of the original equation. When 
one considers the continuum Edwards- Wilkinson equa- 
tion M„ — v{2TTn/NY = vk^. The solution to the set of 
equations for c„ reads 

Cn{t) = c„(0)e-''"* + f dt' e-''"(*-*')??„(0. (9) 



From now on we shall consider the flat initial condition 
Cn(0) = 0, the generalization to arbitrary initial condi- 
tions being straightforward. 

The global roughness of the system is defined as 



N-l 



u{t) 



J=0 



(10) 



which in terms of the Fourier modes can be written as 

^'W = ^ E (cnWc™W)e'('="+'=")^ (11) 



N 

Finally, giving that 

(c„(t)c™(t)) = ^ Az!!^ [1 _ e-K+-)*j , (12) 

N l'n+ I'm L J 

one arrives at the result for the global roughness of the 
system. 



w^{t) = ^y ^ — . 

n—1 



(13) 



which includes the time dependence, the finite size N , 
and the discreteness of the system through f„ given by 
Eq. ©. 

The complete solution of the roughness, Eq. ((T^ . has 
three time regimes. In the context of surface growth, 
where the Edwards- Wilkinson equation was originally 
proposed [3, the three regimes correspond to random 
deposition of particles on a substrate, correlated growth 
of the surface, and roughness saturation due to boundary 
conditions. The different regimes are described by 



2Dt for t < ti 
W\t) = I D.^, for ti «t <Ci2 



% for t2 « t 



with the crossover times given by 

1 



i2 - 



2ttv' 
288v' 



(14) 

(15a) 
(15b) 



The roughness associated to this crossover times are con- 
sequently given by 
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(16a) 
(16b) 



Figure depicts the different discussed regimes for the 
roughness obtained with the DEWE, Eq. ([T^ . The solid 
line represents the roughness of a system with size N = 
128, elastic constant v = 1 and noise intensity D ~ 1/2. 
As can be observed, the three characteristic regimes are 
present, with the asymptotic forms indicated with dashed 
lines. Also shown in the figure are the two crossover 
points yi = {ti.WD and 2/2 = (^2, 1^2^)- 
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B. SFD effective elasticity 

In order to use the DEWE to describe the SFD prob- 
lem, we have to relate its parameters with those of the 
SFD model presented in Sec.|lll Since the noise intensity 
should describe the diffusion of a single particle, it is di- 
rectly mapped to the single particle diffusion coefficient, 
thus setting D = T. Regarding the effective value for 
the elasticity it can in principle depend on both, the 
mean separation between particles L, and the diffusion 
coefficient D. One can argue that v should increase with 
decreasing L, since the particles tend to see each other 
more frequently, thus further restricting the fluctuation 
of its neighbors. It can also be argued that should be 
proportional to the diffusion coefficient, i.e. for larger 
values of D the system becomes more rigid. This is also 
akin to the compressibility of a simple gas, which is pro- 
portional to the temperature of the system. As we will 
show below in detail, the value 



D 



l^SFD — 



(17) 



gives a very good description of the SFD problem in 
terms of the DEWE. Note that in the last equation we 
have used the actual distance between two consecutive 
particles L — 1 (see Fig. [T]), which takes into account the 
size of the particle, that in the present case is equal to 
one. This is important for small separations, but for the 
case L ^ 1, one has that i'sfd ~ 1/L'^. 

One can now use the general results given in the last 
subsection taking D = T and v — j^sfd- The next point 
to be considered is that the flat initial condition in the 
DEWE, Uj{0) = 0, corresponds to a system of particles 
which are initially located in an ordered structure, i.e. 
rj{0) = jL. For future reference in the present work, let 
us recall that using i^sFD, and from Eq. (ITil) . the three 
regimes of the roughness for the SFD problem read now 



( 2rt for t < h 

21^* for ti < t < t2 



(L-1) 

with the crossover parameters given by 



(18) 



27rr 



t2 = 



ttN^{L-1) 
288r 



w. 



2_ N{L-lf 
12 



(19a) 



(19b) 



(19c) 



(19d) 



It follows that for the SFD problem the three rough- 
ness regimes should be interpreted as follow. For t <C ii 
single particles move independently, and thus this regime 
corresponds to normal diffusion. The first crossover time 
marks the time needed for a particle to diffuse over a dis- 
tance of the order of the mean particle separation. Then, 
for times larger than ti particles start colliding with each 
other and the diffusion of the particles becomes corre- 
lated. Therefore, the second regime, ti ^ t ^ t2, cor- 
responds to the subdiffusive motion of particles. In this 
regime, there exists a growing correlation length which 
characterizes the number of particles whose motion is cor- 
related at time t. The crossover time t2 indicates when 
this correlation length reaches the size of the system and 
all particles "know" about the presence of all the other 
particles. Therefore, for t ^ t2 the size of the system 
restricts the growth of the roughness of the system, thus 
corresponding to the saturation regime. 

Figure [3] shows numerical results for the dynamic 
global roughness of the SFD model together with the 
analytical roughness obtained with the DEWE using the 
effective elasticity i'sfd- Each panel corresponds to a 
different number of particles N, and for each N different 
L values are shown. In Fig. [3l symbols correspond to 
numerical simulations of the interacting particle system 
described in Sec.|lll while the continuous lines correspond 
to Eq. (fT5)l with i/ = i'sfd and D — T. As can be ob- 
served, the proposed form for the elasticity does describe 
the numerical data accurately, specially for large values 
of L. 



C. Rougfiness scaling function 

All the different curves shown in Fig. [3j also described 
through the DEWE, can be cast into a single universal 
form depicting the three regimes. When analyzing sur- 
face growth problems, the roughness is commonly scaled 
in terms of the Family- Vicsek form [l^ [l^. However, 
this very usefuU scaling form is well suited for roughness 
presenting two characteristic regimes, as the correlated 
growth and saturation regimes. Here, we are dealing with 
three regimes and thus the Family- Vicsek scaling form 
can not collapse all the data simultaneously. The prob- 
lem of properly scaling a dynamic roughness with three 
regimes has been recently addressed by Chou and Pleim- 
ling (20i] . Based on a phenomenological construction, it 
has been shown that, once the two crossover points for 
each curve are estimated, the scaling can be described 
through 



log (Sf) 

Mm 



= F 



log 



(* 



lot 



(20) 



where F(x) is a given scaling function containing the 
three characteristic regimes. In order to use this scaling 
form the two crossover points must be estimated with 



5 




FIG. 3: (Color online) Roughness of the SFD problem. Points correspond to numerical simulation while continuous lines 
correspond to the DEWE. DifTerent panels are for (a) N = 32, (b) iV = 64, (c) N = 128, and (d) iV = 256. In all panels, the 
different values of L are indicated. 
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FIG. 4: (Color online) Scaling of the dynamic global rough- 
ness shown in Fig. O for different values of the SFD parame- 
ters. 



all the prefactors. It is not enough to only know just 
its functional dependence on the key parameters of the 
problem. 

Figure 2] shows how all the data contained in Fig. [3] 
can be collapsed into a single universal function accord- 
ing to the proposed form, Eq. ((20)) . The data have been 
collapsed using the theoretical values for the crossover 
points, Eqs. p^ . and not the numerical estimates of 
the crossover points. This stresses the relevance of the 
DEWE solution for the SFD problem. 



IV. ADDITIONAL INTERACTION TERMS 

In this section we analyze the effects of an additional 
interaction between particles. We shall consider two dif- 
ferent energetic contributions giving rise to repulsive in- 
teractions, namely a purely harmonic term of intensity 
K and a term including a repulsive interaction decaying 
with the inverse of the squared distance with intensity A. 
It is important to understand the effects of an harmonic 
term since, in principle, any interaction term can be re- 
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FIG. 5: (Color online) Behavior of the global roughness for 
the case of an additional harmonic interaction term and for 
N = 128. The mean separation between particles is L = 32 in 
panel (a) and L = 16 in panel (b), and the different K values 
are indicated. 



FIG. 6: (Color online) Behavior of the roughness for the case 
of an additional interaction term proportional to 1 /r^ and for 
TV = 128. The mean separation is fixed at L = 16 in panel (a) 
and L = 8 in panel (b), with the different A values indicated 
in the keys. 



duced to the harmonic case in a strong interaction limit. 
On the other hand, the other interaction term is com- 
monly used in the vecinal surface system [15] and will be 
analyzed as a typical non-trivial example. 



changing with the additional interaction terms affects all 
the crossover values, as we shal discuss in the next sub- 
section. 



A. Roughness 

Let us first show the general trend of the roughness 
with different interaction intensities. Figure [S] shows the 
behavior of the roughness in the harmonic case and for a 
fixed system size N = 128. The intensity of the harmonic 
term varies in the range 0.001 < K < 0.016 while the 
mean separation between particles is fixed at L = 32 in 
Fig. [Ha) and L = 16 in Fig. [S^b). 

As can be observed, the three regimes are still present. 
The saturation value of the roughness Wf decreases with 
increasing K. In addition, for a fixed K, W2 grows with 
the mean separation L. In fact, this behavior can be ra- 
tionalized by considering that when increasing the inten- 
sity K, the effective mean separation between particles 
decreases since the harmonic interaction generates an ef- 
fective collision before two consecutive particles actually 
touch each other. Therefore, this effective decrease of L 
makes the saturation value of the roughness to decrease. 

Figure [6] shows the behavior of the roughness for an 
interaction term proportional to 1/r^. The number of 
particles is = 128 and the mean separation is fixed 
at L = 16 in Fig. EJa) and L = 8 in Fig. Eljb). The 
intensity of the interaction takes values in the range 
1 < A < 16. The general trend is the same as in the 
harmonic case. For larger intensities, the saturation val- 
ues decrease. This can also be rationalized in terms of 
a decreasing effective separation between particles. The 
fact that the separation between particles is effectively 



B. Crossover points and scaling functions 

We analyze here how the crossover points yi = 
{ti.WD and ?/2 = (i2,W^2^) depend on the interaction 
intensities K and A. For clarity, let us call j/J and j/2 the 
values these parameters take in the pure entropic case, 
given by Eqs. (ITO)) . 

For an harmonic interaction between particles it is 
clear that when K is very large particles will oscillate 
around their initial positions. Then, since two parti- 
cles would hardly be in neighboring sites the hard core 
repulsion can be neglected in the large K limit. The 
system will then be very well described by the DEWE 
with the elastic constant given directly hy v = K. The 
crossover parameters j/i and y2 will take the values given 
by Eqs. ((TS]) and ^ with v = K and D ^T. 

Therefore, the crossover parameters change from the 
K = {) value given by Eqs. to the large K values 
given by Eqs. (TTS]) and ([T5| with v = K . One can expect 
the crossover between these values to be given in terms of 
the quotient K/vqyd^ since this is the effective variable 
which is changing. We then expect that the scaled values 
Vi/Vi be only a function of KL^/D; and the same for 
^2/2/2- We can thus write that 

y--^^-f[-ir)^ (21) 
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with m = 1,2, and the scaling function behaving as 

1 



X 



-1 



for a; ^ 1, 
for X 1- 



(22) 



The x~^ behavior at large K results directly from 
Eqs. and P^ . Furthermore, since the dependence 
of all the crossover parameters with v is the same, we 
expect the same scaling function f{x) for all of them. 

From the data in Fig. [5] we have estimated the 
two crossover points for all curves, which are shown 
in Fig. ^a) according to the proposed scaling form, 
Eq. (j21D . In this figure, close (open) symbols correspond 
to L = 32 (i = 16). The crossover parameters are: ti 
(squares), t2 (circles), (triangles up), and iy| (trian- 
gles down). As a guide to the eye we have also plotted 
in the figure, using a dashed line, a simple function that 
behaves as the scaling function Eg.ip^. This function 
is given by f{x) — (1 + a;)^^. Although the error bar 
of the estimated values of the crossover points is of the 
order of the scatter of all the points in the figure, one can 
fairly observe that the data follow the general trends of 
the proposed scaling form. 

In the case of the l/r^ potential, it is convenient to 
consider the large A limit performing a small u expan- 
sion. To this end, let us consider a simplified situation 
where one particle can diffuse between its two neighbor- 
ing particles which are at fixed positions, and ~L 
respectively. Then, if A is large, the particle will oscil- 
late around w = and the energy can be approximated 
by 



VIA. 
~L4 



— bu 
2 



(23) 



Therefore, in an harmonic approximation, one can con- 
sider the effective elastic constant to be va — 2AA/L'^. 

Since in the large A limit the system can be consid- 
ered within an harmonic approximation, the arguments 
we followed for the harmonic potential can be straightfor- 
wardly used here. In this case, the parameters ym change 
from the value to the values given by Eqs. ([T5|) and 
([T5| but this time with i'sfd = i^a- It follows that 



Vn 



f24A\ 



with the scaling function behaving as 



1 for a; ^ 1, 
^ for x ^ 1. 



X 



(24) 



(25) 



It is also tempting to see if the two functions / and g are 
approximately the same, although this is not obvious a 
priori. 

From the data in Fig. [S] we have estimated the val- 
ues of the crossover points, as shown in the scaled form 
Eq. ([Ml) in Fig.[7Jb). Closed (open) symbols correspond 
to L = 16 (L = 8) and the crossover parameters are plot- 
ted with the same symbols as in Fig. ^a). The dashed 
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FIG. 7: (Color online) Scaling of the crossover points for the 
cases with additional interaction terms, (a) The harmonic 
case, see Eq. (|21|l . with L = 32 (closed symbols) and L = 16 
(open symbols), (b) The 1/r^ case, see Eq. (|24p . with L = 16 
(closed symbols) and L = 8 (open symbols). In both panels 
the crossover parameters are: ti (squares), t2 (circles), Wi 
(triangles up), and W| (triangles down). 



line corresponds to g{x) = (1 + x)~^ and is plotted as a 
guide to the eye. As compared with the harmonic case, 
these data points seem to be closer to the saturation of 
g{x) at very small values. In order to properly reach 
the large A regime one should go to larger system sizes. 
Besides, it should be mentioned here that although we 
expect the large A limit to be described with the scal- 
ing form Eq. ([24]), the small A limit will not necessarily 
be described with this scaling form since the small Su 
expansion is not valid in this limit. 

Finally, since we have already estimated the values of 
the crossover points, we can also present the roughness 
scaling covering the three characteristic regimes by using 
the form given by Eq. ([^0)1 . Figure [5] shows the roughness 
scaling for (a) the harmonic case and (b) the 1/r^ case. 
As can be observed, the data collapse in terms of the 
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FIG. 8: (Color online) Roughness scaling for the (a) harmonic 
and (b) interaction terms. 

estimated crossover points is very good. 

V. UNIVERSAL ROUGHNESS DISTRIBUTION 

As a final characterization of the roughness behavior 
for the SFD problem, we shall discuss here the form of 
the roughness distribution. As has been already shown, 
the roughness distribution can be used as a test of the 
surface growth universality class [22l - [25l |. For instance, 
let be a given instantaneous value of the roughness at 
saturation, i.e. for t ^ t2. The roughness distribution 
P{'uP') can then be written in the scaling form 

- (^) ' (2^) 

where is the average value of the global roughness 
at saturation and thus corresponds to the mean value 
of P{w^). The function $(2;) is a given scaling func- 
tion characteristic of the universality class. The form 




FIG. 9: (Color online) Universal scaling of the roughness dis- 
tribution for the pure hard-core case (A'^ = 64 and L = 8), 
the harmonic potential (A'' = 64, L = 32 and K = 0.008), 
and the potential (iV = 64, L = 8 and A = 6). The 

continuous line corresponds to the analytical solution of the 
Edwards- Wilkinson equation, from Ref. [2l| . 

of <i>(x) for the Edwards- Wilkinson universality class has 
been analytically calculated in Ref. [2l| in the saturation 
regime and in Ref. p6l | in the correlated growing regime. 
In Fig. |9] we show the scaled roughness distribution for 
three data sets corresponding to the purely hard-core, 
harmonic, and 1/r^ repulsion cases. The analytical form 
of $(a;) for the Edwards- Wilkinson case is also shown. 
As can be observed all three cases are very well collapsed 
into the same scaling function according to Eq. (|26l) . 



VI. CONCLUDING REMARKS 

Summarizing, we have shown that the DEWE with an 
effective elasticity i'sfd = D/ {L — 1)^ describes the SFD 
process with hard-core repulsion. We have given a de- 
scription of the process in terms of the roughness, which 
is a global measure of particle fluctuations. With the 
characteristic values of the roughness at the two crossover 
times ti and t2, we were able to scale the roughness in 
terms of a scaling function showing the three regimes. 
Although the phenomenological scaling form works very 
well, we note that it is not clear how it is present in the 
exact solution for the DEWE. 

We have as well explored the effects of different addi- 
tional repulsive interaction potentials between the parti- 
cles. For large intensity values, the dynamics of both the 
harmonic and potentials cases can be described in 
terms of the DEWE with an effective elasticity K and 
VA respectively. We have shown how the values of the 
crossover points scale with the effective elasticity and we 
have properly scaled the global roughness. 

Finally, we have shown that for all cases the roughness 
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distribution in the saturation regime can be recast into 
the universal form of the Edwards- Wilkinson equation. 
This indicates that all the relevant parameters of the SFD 
process are taken into account by the mean value of the 
roughness. 

We have offered here a simple description of the SFD 
process using the roughness of the associated effective 
DEWE, providing an alternative tool to describe the SFD 
problem, which can be complementary to other analyt- 
ical approaches like the ones presented in Refs. 
Furthermore, the use of the effective DEWE opens up 
the possibility of exploring the problem from a new per- 
spective. For instance, related equations also used in sur- 
face growth phenomena, like the Kardar-Parisi-Zhang or 
the MuUins- Herring equations ^T§\, could be potentially 
used to describe processes when the non-passing rule is 
slightly relaxed. In this cases an effective mass transport 
could be considered using these equations. It can be also 
important to consider the effective Edwards- Wilkinson 
equation with quenched disorder to describe SFD pro- 
cesses in a disordered substrate as the one considered in 
Ref. [13! ■ Thus, we expect that the potential relation be- 
tween SFD processes and surface growth equations can 
be further exploited. 

Regarding the experimental relevance of our results. 
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